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Within Tsallis statistics, a picture is elaborated to address self-similar time series as a thermo- 
dynamic system. Thormodynamic-typc characteristics relevant to temperature, pressure, entropy, 
internal and free energies are introduced and tested. Predictability conditions of time scries analysis 
are discussed in details on the basis of Van dcr Waals model. Maximal magnitude for time interval 
and minimal resolution scale of the value under consideration are found and analyzed in details. 
Time series statistics is shown to be governed by effective temperature being exponential measure 
of the fractal dimensionality of a phase space related to the time series. 

PACS numbers: 05.20.-y, 05.90.+m, 05.45.Tp, 05.40.Fb 



I. INTRODUCTION 

Time series analysis allows to elaborate and verify 
macroscopic models of complex systems evolution on the 
basis of data analysis [1]. This analysis is known to be 
focused on numerical calculations of correlation sum for 
delay vectors that allow to find principle characteristics 
of the time series [2]. Being traditionally a branch of the 
theory of statistics, time series analysis is based on the 
class of models of harmonic oscillator which are related 
to the simplest case of the Gaussian random process [3] . 
But well known real time series is relevant rather to the 
Levy stable processes, than the Gaussian ones being very 
special case [4]. Because the formers are invariant with 
respect to dilatation transformation [5], the problem is 
reduced to consideration of self-similar stochastic pro- 
cesses. 

The simplest characteristic of time series is known 
to be the Lyapunov exponent which maximum positive 
value yields predictability domain in the system behavior. 
A range of complexity in such behavior is determined by 
the Kolmogorov-Sinai entropy that equals to sum over 
positive magnitudes of whole set of the Lyapunov expo- 
nents and may be reduced to the usual Shannon value 
in information theory. With passage to nonlinear system 
the probability Pi of i-th scenario of the system behav- 
ior is transformed into power function p^ determined by 
index q < 1, so that the Kolmogorov-Sinai entropy is 
replaced by the Renyi one 



Respectively, governing equation dp/de = —(3p'^, (3 = 
const > describes probability variation with energy e = 
ej to arrive at the Tsallis distribution [6] 

PCX [f - (1 - g)/3e]T^ . (2) 

In the limit q ^ 1, this distribution takes usual Boltz- 
mannian form p oc exp(— /?e) falling down exponentially 
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fast in contrary to the power asymptotic of the general- 
ized exponent (2). Physically, such a behavior is caused 
by self-similarity of the system related to the Tsallis 
statistics [7] — [10]. As a result, the problem appears 
to study self-similar time series that present processes 
corresponded to power-law distribution type of Zipf 
Mandelbrot ones [11]. This work is devoted to analytical 
consideration of such type time series as a thermody- 
namic system. 

The paper is organized as follows. Section II is de- 
voted to elaboration of a model which pcirniits to ad- 
dress a self-similar time series as a nonextensive ther- 
modynamic system. Section III is based on calculations 
of both entropy and internal energy of the time scries. 
As a result, thermodynamic-type characteristics of the 
time series such as temperature and entropy, volume and 
pressure, internal and free energies are introduced. Their 
testing for the model of ideal gas is shown to be basis 
for statistics of self-similar time series. Section IV con- 
tains calculations of corrections to the ideal gas approach 
when external field and particle interaction are switched 
on. On the basis of Van dcr Waals model, we obtain 
the expressions for the specific heat and susceptibility as 
functions of the temperature. Section V is devoted to dis- 
cussion of the physical meaning of the results obtained 
within the framework where the predictability of behav- 
ior of time series is mimicked as stability conditions of 
a non-ideal gas. We find maximal magnitudes for time 
interval and minimal resolution scale of the value un- 
der consideration. Concluding Section VI shows that a 
temperature governing time series statistics is exponen- 
tial measure of a self-similarity index related to fractal 
dimensionality of the phase space. Finally, several equal- 
ities needed in quoting are placed in Appendix. 



II. STATEMENT OF TIME SERIES AS 
NONEXTENSIVE THERMODYNAMIC SYSTEM 

Let us consider d-dimensional time series x(i„) related 
to the set {x„} of consequent values x„ = x(t„) of prin- 
ciple variable x(f) taken at discrete time instants t„ = nr 
that we obtain as result of dividing a whole time series 



2 



length T = Nt by N equal intervals r. It is obviously to 
be relevant to the time series x(<:„), the set {x„} should 
be supplemented by conjugated set {v„} of velocities, 
which show rates of x„-variation with the time jump- 
ing. In the simplest case of Markovian consequence, one 
has v„ = (x„ — x„_i)/r. For more complicated series 
with m-step memory, the velocity magnitude is defined 
as follows: 



Vn = 



(3) 



5i{m,n) 



^{n—m)+i ^(n— m) + (i— 1) 



The paradigm of our approach is to address the time 

scries as a physical system defined by an effective Hamil- 
tonian Ti. = H{x„,v„} on whose basis statistical charac- 
teristics of this series could be found. If one proposes that 
series terms x„ related to different n are not connected, 
the effective Hamiltonian is additive: 



N 



^ — ^ ^ ^n; — ^(Xri, V^^) 



(4) 



Physically, this means that the series under consideration 
is relevant to an ideal gas comprising of A'' identical parti- 
cles with energies e„ . Further, wc suppose different terms 
of time series to be statistically identical, so that effec- 
tive particle energy does not depend on coordinate x„: 
e(x„,v„) e(v„). Moreover, since this energy does not 
vary with inversion of the coordinate jumps x„ — x„_i, 
the function £(v„) should be even. We use the simplest 
square form 



— V 
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(5) 



which is reduced to the usual kinetic energy for a particle 
with mass 1. With switching on an external force F = 
const, particle energy (5) becomes as follows: 



1 



Fx„. 



(6) 



To study a behavior of the time series as a whole one 
needs to fulfil summation over a set of states given by 
manifold {x„,v„} that is relevant to the system phase 
space. In so doing, it is convenient to pass to related 
integrations as following: 



E 



TV 

n 



dx„dv„ 



N 

Y[ JJdyndun. (8) 



{x„,v„} n=l n=l ' 

Here, the factorial takes into account statistical identity 
of the time series terms, A is effective Planck constant 
that determines a minimal volume of the phase space per 
a particle related to a term. The inverted factor 





-dN 




-N 






_N{TAy_ 





Af=N\ 



is caused by change of variables 



Yn 



X 



(9) 



(10) 



rescaled with respect to macroscopic length X being cho- 
sen to guarantee the conditions 



dy„ = 1, / du„ = 1. 



(11) 



According to Introduction self-similarity condition 
forces to use a statistics type of given by the Tsallis dis- 
tribution (2). However, the latter is appeared to be in- 
convenient because it is normalized with nonconventional 
condition 



(12) 



and does not take into account the definition of the in- 
ternal energy 



i \ '9 



(13) 



Finally, when time series has a microscopic memory, di- 
mension of the delay vectors Si{m, n) in the definition (3) 
needs taking m > 1. Moreover, if time series terms x„, 
x„ with m ^ n are clustered, the Hamiltonian becomes 
relevant to a non-ideal gas with interaction Wmn) m ^ n: 



n 



1 ^ 1 



To avoid these difficulties an escort distribution was in- 
troduced [12] 



PL 



(7) In explicit form it reads as follows: 



(14) 



^g{yn,U„} 



(1-9) 



K{y„,u„}-g 



at (l-g)^ 
otherwise. 



.}-E 



< 1, 



(15) 
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Here, the partition function is defined by condition 



Z = N- 



'n//[ 



i-(i-g) 



1-9 



dy„du„, 



(16) 



where < g < 1 is a parameter of nonextensivity, Tg is 
energy scale. Internal energy E is determined by equality 

E = M-^Y[ jj n{yn, u„}P,{y„, u„}dy„du„, (17) 

and normalization parameter (1)^ is expressed by the 
partition function (16) in accordance with Eq.(79). 

To check the statistical scheme proposed let us address 
firstly trivial case of time series x„ = const. Here, the 
particle energy e is a constant as well, so that the Hamil- 
tonian is 7i = Ne. The partition function Z = Af~^ and 
the normalization parameter (1)^ = /^-(^-i) are given 
by inverted normalization factor (9), whereas the inter- 
nal energy E = Ne is reduced to the Hamiltonian. Then, 
the entropy H = —\nj\f obtained according to definition 
(77) given in Appendix is reduced to zero if only the nor- 
malization factor takes the value J\f = 1. As a result, we 
find effective Planck constant: 

--(^)'?- 

Our future consideration is stated on the assumption 
that the volume V = X'^ of d-dimensional domain of 
the coordinate x„ variation in dependence of the particle 
number N is governed by Levy-type law 

X'^ = x'^N-^. (19) 



Here, a; is a microscopic constant and a dynamic expo- 
nent z is reduced to the fractal dimensionality D of self- 
similar phase space [8] 

z = D. (20) 

Then, one obtains the following scaling relation for the 
phase space volume per a term of the time series: 




In the case of Gaussian scattering, when D = 2, the 
minimal volume A'^ of the phase space does not depend 
on number N of the time series terms. Such a condition 
approves our choice of the relation (19) for the whole 
volume V = X'^ as function of the number N of time 
series terms. 



III. NONEXTENSIVE THERMODYNAMICS OF 
TIME SERIES AS AN IDEAL GAS 

Calculations of main thermodynamic quantities of 
nonextensive ideal gas arrives at the following expres- 
sions for the partition function (16) and the internal en- 
ergy (17) (see [12] — [14]) 



^^7(g) 



1-3 



l + (l-g)- 



dN 



l + (l-9) 



E 



(22) 



E = 



2 N\ 



1 + (1 - g) 



i + (i-g) 



E 



(23) 



Here, rf-dimensional gas in volume V 
and the notations are introduced 



is addressed 



to be determined by the Euler F-function. Combination 

of equalities (22), (23) with relation (79) yields explicit 
expression for the normalization parameter: 



e = 



27rr, 
A2 



7(9) = 



(^) 



(24) 



m 



(l + a) 



(25) 
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a=-{l-q)dN. 

In the limits 

1 - g < 2/d, N > 1 

when 

7(5) - [e{l-q)]^ (l + a)-T^?, 
one obtains for the entropy (77): 



Na 



2(1 - a) 



In 



J2+dnd 



(26) 

(27) 
(28) 

(29) 



With accounting scaling relation (19), this expression 
takes the usual form 



G 



(27rer,)^ (-J (30) 



if the dynamic exponent is determined as 
z = D 



1 

1-0/ 



(31) 



Respectively, the internal energy (23) and the normaliza- 
tion parameter (25) read: 



a^^. (32) 



The physical temperature is defined as follows [13] 

W,7^«= f§) ' Ts (33) 



where the last equality takes into account the second of 
relations (32). This definition guarantees the equiparti- 
tion law 



E = CT, C = cN, 



d 



where the quantity 



dE 
df 



(34) 



(35) 



is the specific heat. It is easily to convince that equations 
(30) — (33) arrive at standard thermodynamic relation 



dH _ 1 
'dE^ T' 



(36) 



Above used treatment is addressed to a fixed value of 
the internal energy E [14]. In alternative case when the 
principle state parameter is the temperature T, we should 
pass to the conjugate formalism [12]. Here, standard 
definition 



F=E-TH 



(37) 



F = -CT\n 



of the free energy arrives at the dependence 

T 

Then, the thermodynamic identity 
dF 



dT 



= -H 



yields the relation 



H = C\n 



(38) 



(39) 



(40) 



that plays a role of the heat equation of states. It arrives 
at the usual definition of the specific heat (cf. Eq. (35)) 

C = T?^. (41) 
Let us introduce now specific entropy per unit time 



dN 



to be determined by a minimal entropy 

Hi = {D-l) lnv/27reT« -ln(^) 
and a redundancy 



r = — ln(eA''), 
dr 

Dependencies on the scale x 
D- 1 



D 



a 



iJ 



(42) 



(43) 



(44) 



h{x) = const 



In 



const (45) 



T \T' 

notice that the system behaves in a stochastic manner 

[15]. 

Effective pressure is defined as 

p^-r- (46) 

to measure specific entropy variation with respect to the 
time series scale. Then, we arrive at a mechanic-type 
equation of states 

px^ D-1 (47) 

being additional to the relation (40) of entropy to temper- 
ature. According to Eq. (47), definition of the pressure 
coefficient 

shows that it is fixed by the dynamic exponent (31) being 
the fractal dimensionality of self-similar phase space: 

K = D-1. (49) 

It is worthwhile to note that effective pressure (46) is 
introduced as derivative of the specific entropy h with 
respect to the microscopic scale x. This means a micro- 
scopic nature of so defined pressure which arrives at a 
susceptibility of the time series with respect to a choice 
of the microscopic scale x inverted (obviously, this sus- 
ceptibility is inversion of the pressure coefficient (48)). 
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IV. CORRECTIONS TO THE IDEAL GAS 
APPROACH 



the conditions (27). According to the definition (77), the 
factor (50) yields the entropy addition 



Let us focus now on effect of external field and particle 
interaction which switching on is expressed by equalities 
(6), (7). An external force F = const causes the second 
term in Hamiltonian (6) to arrive at the factor 



sinh(f^) 



FX 
2T 



N 



(50) 



in partition function (22). Here, Zd is a factor depended 
on the dimensionality d only and we put q ^ 1 due to 



i?R.Tt = Na In 



FX 
2T 



(51) 



where we suppress unessential term. With increas- 
ing homogeneous external field, the entropy H^xt grows 
quadratically at FX < T and linearly at FX > T. 

According to [16], cluster expansion of the particle in- 
teraction in Hamiltonian (7) results in additional factor 
in partition function (22): 



2V 



oc 



r(rf/2)' rf^ 



(52) 



where e is effective radius of the particle core. Relevant 
entropy addition 



Hint — — 



2V 



(53) 



monotonously decreases with growth of the particle vol- 
ume V. Increasing temperature T causes the entropy de- 
crease for attractive interaction w < and its increase in 
the case of repelling one [w > 0) . 

Entropy additions (51), (53) arrive at total value of the 
specific heat (41) in the following form; 



C=^N-aN 
2 



FX , fFX\ 

coth - 1 

2T \2T J 



(54) 

2VT ' ^ ' 



where formulae (34) take into account. On the other 
hand, making use of equalities (42), (46), (51), (53) ar- 
rives at the pressure addition 6p determined by the equal- 
ity 



a 


FX , / 






^^coth 1 




d 


2T \ 





N 



(^+|).(55) 



Then, with accounting (49) the pressure coefficient (48) 
takes the form 



K — -I- ''^ < 1 

1 — a d 



FX 
2T 



cosh (^) 



(1 + d)an {v + I) 



(56) 



where effective density 



is introduced. In the limiting case of low temperatures T <C FX/2, we obtain: 

^ d aN , , 
C^-N- — {FX- nw) , 



(57) 



(58) 



a a 
1 — a d 



exp 



FX 



(1 + d)an {v + I) 



(59) 
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Respectively, in opposite case T FX/2, one has: 



2 2T 



{Fxy 

6T 



(60) 



K ~ 



a a 
1 — a d 



1 [Fxy 



r 



d)an (^v 



(61) 



Thus, in the limits of both low and high temperatures, 

the influence of external field F reduces to hyperbolically 
decreasing addition {FX/T)", n = 1,2 to the specific 
heat C; accordingly, the pressure coefScient k gets re- 
ally the constant a/d. In a similar manner, the particle 
interaction affects hyperbolically, as w/T with intensity 
determined by interaction parameter u>, on the specific 
heat, whereas term being proportional to the particle vol- 
ume V is added to the pressure coefficient. 



V. STABILITY CONDITIONS 

It might seem the analysis of a time series would ap- 
peared to be the subject of more much study if both the 
interval t and the scale x take infinitely decreasing mag- 
nitudes (respectively, the number N tends to infinitely 
large values). However, we will show further that exter- 
nal influence and term chistering arrive at predictabilities 
boundaries in behavior of relevant system that are repre- 
sented as instability of modelling thermodynamic system. 
Our analysis is stated on stability conditions 



C > 0, K>0 



(62) 



for the specific heat C and the pressure coefficient k given 
by equalities (54) and (56). 

In the simplest approach being modelled by the ideal 
gas in Section III, we arrives at the natural restriction 
for the dynamic exponent: 



z = D>l. 



(63) 



It means nonextensive dependence of the system volume 
(19) on the particle munber N. On the other hand, in 
accordance with relation (31) the magnitude of the prin- 
ciple index (26) is limited by condition 



a < 1 



(64) 



which restricts the rate of increasing effective density (57) 
with TV. 

As is seen from expressions (54), (58), (60), the spe- 
cific heat falls down monotonically with temperature de- 
crease. To analyze such a behavior quantitatively let 
us consider the case of low temperatures T ^ FX/2, 
when we can use the simplest of above dependencies (58). 



Then, it is easily to see that at temperature less than 
critical magnitude 



FxN- 



wx 



(65) 



the system becomes nonstable due to external force F 
and interparticle attraction w < 0. If the number N is 
much less than critical value Nr defined as 



1 l±iia W - 

TVc = - — -X 

F 



(i+d) 



(66) 



main contribution gets interparticle attraction w < 0, 
in opposite case N ^ N^. — external force F. With 
passage to more complicated case of high temperatures 
T FX/2, when one follows to use temperature de- 
pendence (60), the physical situation does not change 
qualitatively. 

As show estimations (59), (61), the pressure coefficient 
becomes negative if microscopic scale a;'' is less than a 
critical magnitude 



dil + d) /\w\ 



1 



T 



(67) 



that is determined at attractive nature of interaction only 
{w < 0). If both parameters w and v takes nonzeroth 
magnitudes, the temperature values arc characterized by 
another minimal magnitude Ty = \w\/v, lower which the 
microscopic scale Xc becomes principle. This tempera- 
ture is less than critical value (65) if elementary volume 
takes magnitudes lower than value 



d 



d 



-n 



x'^N-" 



(68) 



where we take into account definition (57). 

In the limits v Q, F ^ Q, boundary magnitudes 
(65), (67) arrive at appearance of the critical value 



(l + rf)(l + rf2) 



4d2 



(l + d)(l + d2)2 



(69) 



Relevant fractal dimensionality of the phase space (31) 
grows monotonically with dimensionality of the time se- 
ries d from the magnitude Dc = l + \/2 ~ 2.414 at = 1 
to Dc ~ (i^ at d 3> 1 (see Figure 1). As a result, the criti- 
cal value (69) arrives at the upper boundary of the fractal 
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dimensionality Dc being more than topologic magnitude 
Dt = 2d always {Dc > Dt). Thus, we can conclude the 
upper boundary of application of our approach is ap- 
peared to be inessential at all. 





123456789 



FIG. 1: a. The dependence of the critical value (69) on the 
time series dimensionality, b. The same for related maximal 

value (31) of the fractal dimensionality of the phase space 
(solid line) and topologic magnitude Dt = 2d (thin line). The 
physical domain is shaded. 



VI. DISCUSSION 

We have addressed above the simplest model which 
have allowed us to examine analytically a self-similar 
time series in standard statistical manner. Character- 
istic peculiarity of related equalities is a scale invariance 
with respect to variation of the nonextensivity parame- 
ter 1 — g which is contained everywhere as combination 
(26) derived to the parameter a that is related to the dy- 
namic exponent z and the fractal fractal dimensionality 
D of the phase space according to Eq. (31). This invari- 
ance is clear to be caused by self-similarity of the system 
under consideration. For a given time series, the value D 
is fixed if it is addressed to a monofractal manifold and 



takes a closed set of magnitudes in the case of self-similar 
system relevant to a multifractal. 

In real time scries the property of self-similarity can 
be broken, so that a dependence on the nonextensivity 
parameter itself could appear to be very weak. However, 
above introduced set of pseudo-thermodynamic charac- 
teristics of time series is kept as applicable and accus- 
tomed thermodynamic relations (34) — (37), (39), (41), 
(46) and (48) can be applied to analysis of arbitrary time 
series. 

Main progress in our consideration is that time series 

statistics is governed completely by the temperature (33). 
With accounting Eq.(30) and rescaling the temperature 
unit Tg into Tsc = {2'jTe)"'Ts we derive to the expression 



T 



(70) 



Being independent of the number of terms TV, the time 
series temperature shows exponential dependence on the 
index (26) located under critical magnitude (69). To es- 
tablish character of the power dependence on the ratio 
of the range X of the principle variable to the time in- 
terval T, it is naturally to choice measure units of the 
temperature in the following manner: 



T 



X 



T 



o" " 
(27r)» 



- - 

27r V T 



.(71) 



Then, the expression (70) for the time series temperature 
takes the simplest form 



T = 



D = 



1 



1 



(72) 



according to which the value T is exponential measure 
of the fractal dimensionality D of the self-similar phase 
space related. 

As has shown consideration in Section V, a time series 
subject external influence and term clustering is limited 
in predictability of behavior of relevant system that is 
appeared as instability of modelling thermodynamic sys- 
tem. On the basis of Van dcr Waals model, we have found 
minimal magnitudes (67), (68) for resolution scale Xc of 
the value under consideration. Respectively, making use 
of the critical temperature (65) and the definition (72) 
shows that a time series is predictable if the microscopic 
time interval t takes magnitudes less than a critical mag- 
nitude Tc- In the case N ^ N,., where the boundary 
magnitude is determined by Eq. (66), we find 



/a, ,\ 2 D -I I d 1 

In opposite case A'' Nc, one obtains 

_ 1 

/ a \ 201 ^ ^ 1 

Tc = (^-Fj e 2 X 2 iV 2-D<l . 



(73) 



(74) 



It is easily to convince that the entropy (40) takes positive 
values within the stability domain t < Tg. 
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To conclude, we point out that one of advantages of 
the approach proposed is a possibihty of its apphcation 
to numerical analysis of real time series. The basis of 
this calculations is stated on expressions (3) — (7), (15) 
- (17), (33), (34), (35), (37), (39), (41), (42), (46), (48) 
and (77). This work is in progress. 
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Appendix: Main statements of the nonextensive 
statistics 

Our approach is stated on using entropy definitions 

H = a\n[expg{Hq)], if, = lng[exp(if/a)]; (75) 
a = ^{l-q)dN 

that are alternations of the usual functions logarithm 
ln(a;) and exponential exp(a;) with corresponding Tsal- 



lis generalizations [6] 

— 1 _j_ 
\ng{x) = , expg{x) = [1 + (1 - q)x] i-^ . (76) 

In explicit form, relations (75) are appeared as 

H = alnZ, Ha = ^^'l" ~ ^ . (77) 
I — q 

Making use of the first equality (75) and the second for- 
mulae (76) and (77) shows that physically defined en- 
tropy H is extensive value being reduced to the Renyi 
definition (1): 

H = aK, = ^ln{l}^. (78) 

Comparison of this result with the first of equalities (77) 
arrives at the relation 

(1), = Z'-'^ (79) 
that ensures the normalization condition [6]. 
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